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Abstract 

We study an integral equation that extends the problem of anti-differentiation. We formulate this equation 
by replacing the classical derivative with a known nonlocal operator similar to those applied in fracture 
mechanics and nonlocal diffusion. We show that this operator converges weakly to the classical derivative 
as a nonlocality parameter vanishes. Using Fourier transforms, we find the general solution to the integral 
equation. We show that the nonlocal antiderivative involves an infinite dimensional set of functions in 
addition to an arbitrary constant. However, these functions converge weakly to zero as the nonlocality 
parameter vanishes. For special types of integral kernels, we show that the nonlocal antiderivative weakly 
converges to its classical counterpart as the nonlocality parameter vanishes. 

Keywords: nonlocal operators, integral equations, Fourier transform, tempered distributions, 
peridynamics, nonlocal diffusion 


1. Introduction 


We consider integral equations for distributions u e on R of the form: 


D a ^u e (t) := —a e * u e (t) = F(t ), t e 


( 1 . 1 ) 


where * denotes the convolution, f * g (t) = J R f(t — s)g(s)ds, a e is an anti-symmetric function on R that 
depends on a positive “nonlocality parameter” e, and D a e is a “nonlocal derivative”. 


Nonlocality describes interactions over distances; in (1.1), it refers to the fact that F{t) is related to u(s) 


via D for values of s far from t, where “far” is quantified by the parameter e. Mathematically speaking, 
D ae is “strongly nonlocal” in the sense of Rogula [T], since the support of D a e u is not contained in that of 
u. In contrast, the derivative u'(t ) does satisfy this property, so it is a “local” operator (it is also “weakly 
nonlocal”; see CD- 

We define D a t in such a way that D a e f(t) —> f'(t) in some sense as e —> 0. In other words, we think of 


(1.1) as a “nonlocal extension” of the first order ordinary differential equation (ODE): 


u' (t) = F(t ), t £ R. 


( 1 . 2 ) 
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We solve (1.11 using the Fourier transform: 


/ OO 

e~ 27Ilit u(t)dt, 

-CO 

/ OO 

e 27rit «u(Ode 

-OO 


(1.3) 


Nonlocal models, though well known since the 1800’s viz. fractional derivatives, have recently found 
successful application in Silling’s [?' theory of peridynamic fractures; for applications of these nonlocal 
operators in nonlocal diffusion and image processing, see Oil!. 

The success of nonlocal models stems from the ease with which they handle singularities. Although clas¬ 
sical models using differential equations cannot well describe discontinuous solutions to these equations, such 
as those occurring in fracture dynamics, the integral operators applied by Silling suffer no such difficulties. 
Using nonlocal operators in these classical models extends the types of physical processes that we can model 
and the types of qualitative behavior that we can describe. 

An important feature of these operators is that they are extensions of classical differential operators. In 
[5], it is shown that such first order nonlocal operators converge strongly in L 2 to classical partial derivatives 
as a nonlocality parameter vanishes. This extensivity property is important since it allows us to preserve 
much of the physical structure that makes up classical differential models. In addition, when the classical 
descriptions are correct, the nonlocal frameworks can recover these results by passing to the “classical limits”. 

The majority of the nonlocal literature has focused on so-called second order models. These nonlocal 
models extend second order partial differential equations and boundary value problems, such as the wave 
equation [6], Laplace’s equation 0 and the heat equation [5] (see also the nonlocal counterpart to the fourth 
order biharmonic equation 0). However, the literature for nonlocal extensions of first order models is more 
sparse. Du et al m studied a nonlocal extension of the nonlinear advection equation. Among other results, 
they showed that inviscid solutions to this nonlinear equation do not blow up in finite time, a stark contrast 
to those of the classical equation. There is also some literature available for first order models using fractional 
derivatives (see m for a large set of examples). 

Given the success of these nonlocal operators in extending second order models, we are interested in 
their application to first order models, specifically first order ordinary differential equations (ODEs). First 
order models are ubiquitous in the natural sciences. One area that is currently under active research is 
the treatment of discontinuous models. These are ODEs that have discontinuous “forcing terms”, and have 
applications in the flow through porous soil [12] . static friction problems [ 133 , an d optimal control with 
discrete feedbacks El. A challenge with classical (i.e. differential) frameworks is to solve these problems 
numerically [IS], since the classical derivatives for solutions to these problems do not exist everywhere. As 
a result, many heuristic and complicated numerical methods are needed to solve these equations computa¬ 
tionally. Although we do not pursue this here, one motivation for studying nonlocal first order models is the 
possibility of replacing classical derivatives in these equations with nonlocal derivatives. It is possible that 
solving discontinuous integral equations, which do not contain classical derivatives, is a more straightforward 
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task than solving discontinuous ODEs. 

The contributions of this paper are as follows. 


• Weak operator convergence: We show that the nonlocal derivative converges weakly to the classical 
derivative. 


• General solution for antiderivative: Using Fourier transforms, we find a general formal expression 
for the nonlocal antiderivative. In general, it is infinite-dimensional, but we identify a class of kernel 
for which it is one-dimensional, analogous to the classical solution. 

• Weak convergence of solutions: Given some hypotheses, we show that the nonlocal antiderivative 
is a tempered distribution and that it converges weakly to the classical antiderivative. We also show 
strong L p convergence under additional hypotheses. 


In § [2j we give an overview of the nonlocal derivative operator and prove that it converges weakly to 
the classical derivative. Section [3] focuses on nonlocal integration. We find the general form of the nonlocal 
antiderivative in § |3.1[ provide explicit examples in § |3.2| and discuss the convergence of the antiderivative 
to the classical solution in § |3.3| After showing in § |3.3.1 that the homogeneous solution converges weakly 
to the constant function, we prove the inhomogeneous term’s convergence in § |3.3.2| 


2. Nonlocal derivative 


The nonlocal derivative operator used here was presented before in [5]. 

Definition 2.1 (Nonlocal derivative). A nonlocal derivative D a e is a convolution operator acting on u : 
M —» R with kernel a e depending on a nonlocality parameter e > 0: 


D a ^u{t) := —a e * u(t) = / a e (s — t)u(s)ds, t £ 

Jr 

The kernel should satisfy the following conditions: 
(i) Scaling and normalization: 

p i-i 


( 2 . 1 ) 


a e (s) = a e a(y/e), a e = e 2 
(ii) Finite dipole moment: 


e / sa(s)ds 

R 


0 < / sa(s)ds 

Jr 

(iii) Anti-symmetry: 
a(s) = sp(|s|). 


=: |a (1 )| < oo, 


( 2 . 2 ) 


(2.3) 


(2.4) 
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Remark 2.1. We will always impose another condition that regulates the behavior of a at infinity. How re¬ 
strictive this condition needs to be depends on the application. Examples: suppa C [— 1,1], or f R |s J a(s)|ds < 
oo for each j = 1 , 2, .... 

Remark 2.2. Du et al [5] used a weaker form of normalization than the explicit scaling in 
convenience of this definition may be worth the cost in generality. 

These conditions characterize D at as a nonlocal extension of d/dt. They ensure that D a e u —> v! in some 
sense as e —> 0. Of course, in which sense this convergence occurs and the space that D a , e u lives in depend 
greatly on the space that u lives in. For example, Du et al [5] proved that, if a has compact support, and u 
is in Wf, then D a ^u is in L 2 and that || D a e u — u , ||/ / 2 —> 0 as e —> 0. But even if it is a distribution (i.e. a 
bounded linear functional on the space Cf°), we can still obtain weak convergence. 


(2.2), but the 


Theorem 2.1 (Weak convergence). Let u be a distribution and a have compact support. Then D a e u —*• v! 
in the Cf° topology. 


Proof. We first show strong convergence for any if in C° 


D, 


t ,eip(t) = / ajs)ip(s + t)ds 

J R 


E (*) 


V>(t) + sip'ft) + -s 2 ip"(6(s, t)) 


By Taylor’s Theorem and (2.2): 
d s 


I loo [ s 2 |a e (s)|ds 
z JR 

= i>'(t) + i^ (2) i e 11 V ,,/ 11 oo ? 


(2.5) 


where ||.||oo denotes the supremum norm, and 0(s,t) £ (t, s + t ) comes from Taylor’s theorem. It is clear 
from this estimate that || D a ^ip — ip '\\oo —t 0 as e —> 0. 

Since ^D at€ ip = D a ^ip ’, and ip' is in Cff 3 , we see that the above result holds for any derivative of ip. 
More precisely, for any n = 0,1,2,..., we have || Da^ip^ — ip( n+1 '> < e|a|( 2 ) ||i/^ n+2 ) ||oo/2|a:(i) | —> 0 as 
e —> 0. 

Since a has compact support, D a ^ip is in C£° for each e, with supp D a ^ip C suppD Qi i if e < 1. Thus, 
D a ,eip —»■ ip' in the (strong) sense of Cf° (for more details on this topology, see e.g. jISlH7j ). 

Let (,) be the standard inner product. Since u is a continuous functional on CfP, we can use the strong 
convergence to conclude that (u,D a , e ip) —» ( u,ip') for any ip in C%°. But we define the distributions D a e u 
and u' by the inner products ( D at u,ip ) := (it, D* t ip) and (u',ip) := —(u,ip’), where D* a e is the adjoint of 
D a ^ e . In fact, as in the classical case, we have D* at = —D a ^ [Sj. This means that ( D a ^ e u,ip) — > (u',ip) for 
every ip, or that D a e u v!. □ 


Remark 2.3. We assumed that a e had compact support in order for (it,a e * ip) to make sense, given that u 
is a general distribution. Different assumptions about the space that u lives in require different restrictions 
on the kernel a. But in each case, the proof of weak convergence proceeds almost identically. 

By imposing stronger hypotheses on it, it is likely that this weak convergence result can be replaced 
by something stronger, as Du et al. [5] showed for Wf functions it. Figure [l] demonstrates convergence for 
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u(t) = \t\ 1 / 2 and a(s) = sl ( _ 1 . 1) (s), where l(^i,i) is the indicator function on (—1,1). The weak convergence 
result in Theorem |2.1| applies, but the convergence appears to be quite strong. Indeed, away from t = 0 
(say, on R\ [—1/2,1/2]), the convergence of the nonlocal derivatives to the classical derivative appears to 
be uniform even in the higher derivatives. The reason for this is that u is a continuous function that is C°° 
almost everywhere, so it is significantly more regular than a distribution. 


£Wu(t) 



e = 0 
e = 0.1 
e = 0.5 
e = 1 


Figure 1: D a ^u(t) for u(t) = \t] 1 / 2 and a(s) = sl(_ 1[1 )(s). 


We note that the anti-symmetry condition (2.4) is critical for interpreting D a e as an extension of the 


classical derivative d/dt, or as a “first order” operator, as evident from its role in (2.5) of canceling the u(t ) 
term. Symmetric kernels a e would instead give something similar to a “second order” operator like d 2 /dt 2 ', 
indeed, convergence to such for a similar integral operator was proven in j!8j . 


3. Nonlocal integration 

In the classical setting, the problem u'(t) = F(t) yields, in general, the solution u(t) = A + f F(t)dt (i.e. 
the “antiderivative”). We investigate this problem using the nonlocal derivative D a e in place of d/dt. More 
specifically, we consider a given function F = F(t) and an unknown distribution u e that satisfy: 

D a ,e'Ue{t) =—a e * u e (t) = F(t), f SR. (3.1) 

A similar type of convolution equation has been considered in harmonic analysis for F = 0 [19], f20l . 
in which case u is said to be a “mean-periodic” function. See also m for the case of nonzero F when 
u(t < 0) = 0. In these studies, the kernel is usually taken to be a measure, so it satisfies a non-negativity 
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requirement. We work with an anti-symmetric kernel, so our work is more analogous to the study of Hilbert 
transforms. The main distinction of our study, however, is the connection of these integral equations with 
the classical differential equation in the limit of e —» 0. 

Our results can also be connected to the distribution of zeros of a e , the Fourier transform of a e . For 
entire analytic functions in the complex plane, there is an extensive theory developed in this area; we cite 
the seminal work of Polya [223 and the review article by Dimitrov and Rusev [23]. We also mention the 
works by Sedletskii, summarized in [M], that expand on Polya’s work and also study the completeness of 
the exponential basis functions arising from these zeros. 


3.1. General solution 


In general, we can find a formal solution to (3.1) by taking its Fourier transform (1.3): 


-a e u e (£) = F(£), £ G R. 


(3.2) 


We can formally solve this equation by division if we make some assumptions about how or if a e vanishes. 


Definition 3.1. We consider those kernels a that satisfy conditions (ii)-(iii) in Definition 2.1 as well as the 
following condition on its behavior at infinity: 


[ |s- 7 a(s)|ds < A a j\ 

Jr 


(3.3) 


for some A a > 0 and each j = 1, 2,.... 


Remark 3.1. Examples of a that satisfy this strengthened condition include rapidly decreasing functions like 
a(s) = sgn(s) exp(— |sp) or those with compact support, such as a(s) = s/?(|s|)l(_i.i)(s). 


If a satisfies conditions (2.3), (2.4), and (3.3), then a is entirely analytic on the real line. This means 


that a e has (at most) countably many real zeros, a finite number of real zeros in each bounded interval, and 
zeros of finite multiplicity. See e.g. [23] for examples of stronger conditions that give more precise control 
over these zeros. 


Given this assumption on a, we find that the general (formal) solution to (3.2) is: 

s(fc) 


«<=(£) = - F (£)/&e(Q + 


E 


A jtk S^(0, 


(3.4) 


£(U 0 l«) m> (Gdl=o 


where g^ is the fcth derivative of g , S a {x) = Sq(x — a) is the Dirac measure centered at a, and each Aj^ is 
an arbitrary constant. The lower sum selects over roots £ j jlf of de(£) with multiplicity k + 1. 

We interpret this expression as a formal infinite series. Note that, by the anti-symmetry condition 


(2.4), we have that d e is also anti-symmetric, so the Ho : o So term always exists. We order the zeros like 
• • • < < £-n+i <■••<£() = 0<---< £„_i < £ n < .... By the anti-symmetry of a e , we have 

€—j,e = — £j,e- 
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Applying the inverse Fourier transform to (3.4) and renaming the arbitrary constants Aj k gives the 
following formal solution to (3.1): 


u e (t) = -T~ 1 [F/a t }(t) + (3.5) 

E^ = 0 l«' m) (O,e)l=0 

Note that the j = k = 0 term corresponds to A 0 0 = constant, which also occurs in the classical antiderivative: 
u(t) = f F(t)dt + A 0 ,o- The first term in (3.5) is analogous to the integral of F, and is exactly such if 
a e = —27r*£. 


3.2. Examples 

We are interested in the nonlocal derivative D a e being an extension of the classical derivative d/dt , so 
we would like to know for which kernels a these operators have similar properties. 

The classical antiderivative u(t) = A + f F(t)dt only involves one arbitrary constant, while the nonlocal 
antiderivative (|3.5|) involves, in general, an infinite number of such constants. We observe that requiring 


|d^(0)| > 0 and |&(£)| > 0 when |£| > 0 results in all of the A^ k in (3.5) vanishing except for A 0j0 . The 


first property holds by (2.3). We can guarantee the second property by requiring that a(s) be a decreasing 
positive function for s > 0; for example, a(s) = sgn(s)/3(|s|) for any /3(|s|) decreasing with |s| satisfies this 
condition; of course, there are many more functions, such as a(s) = sexp(— |s|), that do not need to be 
decreasing to satisfy this condition. For a similar problem on the characterization of the positivity of Fourier 
transforms, see [251. 


Thus, for these types of a e , (3.5) becomes: 


u e (t) = A - F 1 [F/d e \(t), 

where A is an arbitrary constant. This is much more analogous to classical integration. 

As an explicit example, consider a e (s ) = sgn(s) exp(—|s|/e)/2e 2 and F(t) = 1/(1 + t 2 ). The classical 
antiderivative of this function is obviously u(t) = A + arctan(f). On the other hand, for this a e , its nonlocal 
antiderivative is u e (t) = u(t) + 2e 2 t/(l + t 2 ) 2 . We clearly have convergence of u e to u as e —>• 0, so this 
nonlocal antiderivative extends the classical one. 

For another example, consider a(s) = sin(7rs)l^_i.i)(s). Then the Fourier transform of a e is a e (^) = 
*sin(27re£)/[e(4e 2 £ 2 — 1)], which is zero at £ = 0, ±l/e, ±3/2e, ±2/e,.... Indeed, one can directly check with 


integration that — a e * g n = 0 for g n (t) = exp(nirit/e) and n = 0, ±2, ±3,..., which agrees with (3.5). So, 


in addition to annihilating constants, the nonlocal derivative D a t for this kernel actually annihilates an 
infinite dimensional set of functions. This is a stark contrast to the classical derivative. Note, however, that 
as e —>• 0, each nonzero £j jC tends to infinity, so, in a weak sense, the nonlocal antiderivative still recovers the 
classical antiderivative. 
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3.3. Convergence to classical solution 


We demonstrate that (3.5) converges weakly to the classical solution of v! = F, namely u = A + f F. 


There are two types of terms in (3.5) that we need to consider: the homogeneous terms (summation) from 
the equation D a e u e = 0 and the inhomogeneous term that involves F. 

3.3.1. Homogeneous part 

We first show that the zero set of a e spreads away from the origin as e —>• 0. The scaling assumption in 


(2.2) is critical here. 


Lemma 3.1 (Zero Set). Let a e satisfy the conditions of Definitions 2.1 and 3.1 Let denote the jth zero 
of a e as in (3.4). Then, for each j = 1, 2, 3,..., we have —> ±oo as e —> 0. 


Proof. Since f-j. e = — £j, e , it suffices to show that —> oo for each j = 1, 2, 3,.... At a zero £y e , we have: 

/ OO 

sin(2Tr£j je t)a e (t)dt. 

-OO 


Using the scaling in (2.2), we can rewrite this after letting t —> es as: 

/ OO 

sin(27re£y e s)a(s)ds = d(e£y e ) =: 

-OO 

where fj := e£y e is the jth root of a and does not depend on e. Thus, /e —> oo as e —> 0. 


(3.6) 


□ 


We can now show that the homogeneous part of (3.5) converges weakly to the classical solution (i.e. the 


constant function) as e —>■ 0. We first consider finite sums in (3.5). 


Theorem 3.1 (Finite sums). Let u t be the solution of (3.1) given by (3.5) forF = 0, and let a t satisfy the 
conditions of Definitions \2~t\ and 3.1. Suppose that the formal sum in (3.5) is finite. Then u e — Ao,o as 
e 0. 

Proof. Let ip be an arbitrary function in Cf°. For any j : |j| > 1 and any k > 0, the Riemann-Lebesgue 
show that (t k e 2m ^ j ’ et , ip) — > 0 as e —> 0. Since ip is arbitrary, it follows that each 


lemma and Lemma 


3.1 


nonconstant term in the formal sum ( |3.5| ) converges weakly to zero. The result follows since the sum is 
finite. □ 


We can extend this to infinite series if we restrict the formal sum in (3.5) to converge in a suitable 


sense. It is convenient to map the terms of the series to the sequence { U nj£ } and to consider the partial 
sums UN,e = Y2n= i U n ,e- To interchange limit operations, we also make use of the counting measure dp,: 
I fndpn ■= En=l/«- 


Theorem 3.2 (Infinite series). Let u e and a e be as in Theorem 3.1 Suppose that the following hypotheses 
hold: 

(i) For each e > 0, UN,e(t) —> u e (t) a.e. in K. as N —> oo. 


















(ii) For each e, the sequence {wiv,e}“ = i of partial sums is uniformly locally integrable. That is, for every 
n and each £c8 with \E\ < oo, we have \J E UN,e(t)dt\ < U e (\E\), where U e > 0 is locally bounded, and 
Ll e {x) —> 0 as x —> 0. 

(in) For any ip in Cf° , the set of integrals {(U n ^,ip)} e> o is uniformly bounded for almost every n by a 
summable sequence {v n ^} C R+. That is, for each e > 0 and each ip in C^° , we have v n ^ > \(U n>t ,ip)\ for 
every n = 1,2,3,... except for a finite set, and Y^^Li v n,ip < oo. 

Then u e A 0 0 as e —> 0. 


Proof. Fix an arbitrary ip in C£°. Since |supp "01 < 00 and (ii) holds, Vitali’s convergence theorem applies 
to {wAr, e }^ =1 . We compute: 

N r 

lim(tt e ,^>) = lim ( lim un e,ip\ = lim lim N ( U ne ,ip ) = lim / (U n t ,ip)dp n . 

e—>-0 e—>-0 \N—>oo ’ / e—>-0 N—>oo ^' ’ e—>-0 / 


n—1 


By (iii), the dominated convergence theorem applies to the last integral, so by Theorem 3.1 we obtain: 
lim(it e , ip) = / lim {U ne ,ip)dp n = / 0 dp n = 0. 


□ 


Remark 3.2. If we allow the coefficients A :h p to vanish with e, then this weak convergence can be made 
strong (say, in L% c ). 


Remark 3.3. Condition (ii) is satisfied if UN,e —> u t in L) oc . Condition (iii) is satisfied if the zeros of d e have 


a maximum multiplicity of M > 1 and if y~C k \Aj^\ < oo, where the coefficients Aj^ are as in (3.5). Indeed, 
since U ntE (t) = A n t k "for some j n , k n , and A n =: Aj n ^ n , we have |I/ njfr (f)| < |A„||f| M for each e. 
This gives: 


£|(tW)|<£K| [ \t M m\dt<2T»f +1 \\iP\\ 00 J2\An\, 

n n J supp ill n 


where = max( | inf supp ip\,\ sup supp ip\). 


3.3.2. Inhomogeneous part 


We now investigate the convergence of the inhomogeneous term in (3.51 that depends on the “derivative” 


F. 


Definition 3.2. We denote the inhomogeneous part of u e in (3.5) as follows: 


v e (t) = -J 7 1 [F/a e ](t). 


(3.7) 


We first show under what conditions this formal solution v e is well-defined. Then, provided that it is 
well-defined, we show that v e —<■ v := f F, or that it converges to a classical antiderivative in the weak sense 
of tempered distributions. 

A necessary condition for v e to be well-defined is for F to make sense. We recall the definition of a 
tempered distribution. 
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Definition 3.3. Let S be the Schwartz space of infinitely differentiable rapidly decreasing test functions. 
Equivalently, ijj is in S if it is in C°° and if limiq^oo \t\ n ip(t) = 0 for every integer n > 0. The space of 
tempered distributions S’ is the continuous dual to S (i.e. the set of bounded linear functionals on S). 


Definition 3.4. We require that F in (3.1) and (3.5) be a tempered distribution, and we allow it to be a 
generalized function. 


Remark 3.4. Unfortunately, this is a very strong condition. In the classical case, we could find an antideriva¬ 
tive of any locally integrable function F (say, F(t) = e 4 ). A similar problem concerns solutions to An = 0; 
there are many solutions to this equation, such as e Xl COSX 2 , that are not tempered distributions, and Fourier 
transform techniques cannot recover these solutions. 


Remark 3.5. Note that, from this restriction and Definition 


3.1 


makes sense, since we have that 2ni/a e — l/£ is a C°° function, 
makes sense as a tempered distribution. 


the distributional product F/a e always 
Of course, we still need to show that it 


We now wish for the inverse Fourier transform in (3.7) to make sense. There are two problems with this. 
The first problem concerns the possibly infinite number of singularities of 1/oy, and the second problem 
comes from the behavior of 1 /a e at infinity. If each of these behaviors is not sufficiently controlled, then 
F/a t fails to be a tempered distribution, and ( |3.7| ) makes no sense. 

To solve the first problem, we could impose that 1/d have a finite number of singularities. However, for 
simplicity, we will restrict 1/d to be singular only at the origin. The more general case proceeds similarly. 
Note that this means we no longer need the full strength of Definition |3.1[ but we keep it for simplicity. 


Definition 3.5. We require that d(£) = 0 for £ in R. only if £ = 0. A sufficient condition for this (in addition 


to those in Definitions 2.1 and 3.1) is as follows. We require that a(s) = sgn(s)/3(|s|), where /3 > 0 is a 
decreasing function that is strictly decreasing a.e. in (0,a a ) for some a a > 0. 


Remark 3.6. The condition is sufficient by the alternating series test. Indeed, the Fourier integral J 0 °° sin(27r£s)a e (s)ds 
X^o( — l) n Jo 72 " sin(27r£s)a e (s + n/2£)ds > 0 if £ > 0 and if this condition is satisfied. Examples of admis¬ 
sible (3 include (1 — s m )l( 0 ,i)(s) for m > 0, exp(—|s| m ) for m > 0, and |s| _m l( 0! i) for — 2 < m < 0. 

As it turns out, this condition almost solves the second problem as well. We only need two modifications 
to it. We gauge the behavior of a near zero using power functions, and we impose that, near zero, a cannot 
be flatter than a continuous power function. 


Definition 3.6. There are two cases to consider, depending on how singular a is at zero. Let «o+ = 


linisx^o a(s). In all cases, we presume the existence of some endpoint b a : 0 < b a < a a (cf. Definition 3.51, 
some exponent k a that satisfies —2 < k a ^ 0, and some constants K± a > 0. 

(i) If ao+ < oo, then we require that —A'_ jCt s fc “ < a(s) — ao+ < —K+,aS ka for each s in (0,6 Q ). Here, 
0 < k a , and K- a > K+, a - We assume that b a can be made small enough so that we can choose K±, a that 
satisfy K +ta (s + l) fc “ — AT_ iQ ,s fe “ > 0 for each s on (0,1). 
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(ii) If «o+ = oo, then we require that /\_ jC ,s fc “ < a(s) < /v +;Ct s fe “ on (0,b a ). Here, —2 < k a < 0, 
and K-. a < K+ : a- We again assume that, by choosing b a small enough, we can find K±, a such that 
K_ a s ka - K +ta (s + l) fe “ > 0 on (0,1). 


Remark 3.7. “Flat functions” like a(s) = sgn(s)(l — exp(— |s| 2 ))l( 0 ,i)(|s|) are not admissible. 

Remark 3.8. By anti-symmetry, analogous inequalities hold for s < 0. Although power functions are not the 
most general gauge functions (e.g. logarithms are also natural choices), they are quite convenient here due 
to their homogeneity under dilations. 


Remark 3.9. Since a satisfies both the anti-symmetry and positivity conditions of Definitions |2. 1| and 3.5 


we see that it cannot be continuous at zero. It either approaches a constant or blows up as s \ 0. This 
means that ao+ > 0. 


To prove that v f in (3.7) is a tempered distribution, given these restrictions on a, we need the following 
two estimates on a e . They elucidate its behavior for small and large |£|. 


Lemma 3.2 (Near-field). Let a e satisfy the conditions of Definitions 2.1 3.1 and 3.5 Then the following 
estimate holds for each £ in M: 


2tt]£| - C a e 2 1£| 3 < *d e (|£|) < 2 tt|£|, 

where C a = 47r 3 a( 3 )/3a(i), and ay) = f^° s^a(s)ds. 


(3.8) 


Proof. Using the anti-symmetry of a, we can write the Fourier transform (1.3) of a e as a Fourier-sine 
transform: 


/•OO 

—a e (f)/2i= / sin(27r£t)a e (£)dt. 
Jo 


Using the scaling in (2.2), we can rewrite this after letting t —>■ es as: 

1 


~<M£)/ 2 * = 


ea (i) Jo 


sin(27re£s)a(s)ds. 


(3.9) 


Suppose first that £ > 0. Using that sin(27re£s)a(s) < 27re£sa(s) in (3.6) (note the positivity from Definition 


3.5) and that / 0 °° sa(s)ds = a(i)/2 gives: 


*d e (£) < 2 tt£. 


On the other hand, using Taylor’s theorem, letting 6(s) be in (0, oo), and applying the estimate cos(0(s)) < 1 


to (3.6), we obtain: 


fd e (£)/2 = tt £ - 


47r 3 £' J e 
3a 


3,2 r °° 


( 1 ) Jo 

3„, t3,2 , 


cos(0(s))s 3 a(s)ds 


> tt£ - 27T 3 a (3) £ 3 e /3a (1) . 


If £ < 0, then we reverse the above inequalities. This gives (3.8). 


O 
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Lemma 3.3 (Far-field). Let a e satisfy the conditions of Definitions \Kl\ \3.1\ \3.5\ and 3.6 Then, for each 
e > 0 and each £ : |£| > 2eb a (cf. Definition\3.&\), the following estimate holds: 


where k a is as in Definition \3.6\ and C' a > 0. 

Proof. Suppose that f > 0 (by anti-symmetry, we can reverse the resulting inequality when £ < 0). We use 


(3.10) 


the positivity in Definition [3^5] and the scaling in Definition 2T to obtain the following alternating series-type 
estimate: 


nOO 

ia e (£)/ 2 = / sin(27r£s)a e (s)ds 

Jo 


-i/e 


> 


sin(27r£s)a: e (s)ds 


1 


(3.H) 


«(i)e 2 £ 

-JL_ / 

2a(i)£ 2 C Jo 


sin(27rs)a(s/e£)ds 


sin(7rs)[a(s/2e^) — a((s + l)/2e/£)]d s. 


We estimate each term using the inequalities in Definition 3.6 To do this, we suppose that £ > 1/2 eb a (cf. 


Definition 3.6). For each case, (3.11) becomes: 


Case (i): 
ia e (£) > 


2 k c 


a (1) e 


2 +k c 


£ 1+fec 


■ f sm(ns)[K+ a (s + l) ka — 

Jo 


^2 -\-k c 


Case (ii): 
ia-eiO > 


£l+fco' 


1 


2 fea a(i)e 

^OL 


2 +k c 




■ f sin(7rs)[A'_ :Q s fe “ — K +tCl (s + l) fc “]ds 

Jo 


^2 -\-k c 


^1+fco. ‘ 


□ 


Theorem 3.3. Let F satisfy Definition \3.4\ and a e Definitions roifjTifxi and \3.6[ For each e > 0, the 
generalized function l/a. e is a tempered distribution. In addition, v e = —F/a e is a tempered distribution. 

Proof. We let l/a e (£) = —1/27ri£ — /3 e (£)/27r i, where, for |£| sufficiently small, (3 e satisfies the following 
estimate (cf. ( |3.8| )): 

C a <?\£\ 


= |2tt£- m e (g)| < Cqe 2 !^ 2 < 


1^(01 


<*(01 - (27T - C' a e2|e| 2 )' 


(3.12) 


This shows that $ £ is continuous at £ = 0. Since a e has no other zeros (cf. Definition 3.5), we conclude that 
/3 e is a locally bounded function on R. (i.e. it has no singularities). 
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We claim that is a tempered distribution. Since it has no singularities, we need only examine its 

in terms of /3 e (£) = —1/£ + 2n /id e (£), we obtain the 


asymptotic behavior at inhnity. Using Lemma 
following inequality for |£| > 1/2 eb a : 

„ 2lTf 2 + k <X 1 

ACKDIS —ICI-H— 


3.3 


(3.13) 


Since this holds for each e > 0, we see that the function /3 e is a tempered distribution. Since l/£ (interpreted 
in the principal value sense) is a tempered distribution, this means that 1 /a e is also a tempered distribution. 
To show that v e is a tempered distribution, we write the product v e as follows: 

2niv e = —2 tt iF/a e = F/% + Ff%. 

If we interpret l/£ in the principal value sense, then the first term corresponds to the distributional integral 
of F (see e.g. Problem 6.9 in HZ] ). The second term is also a tempered distribution, since (F s 0 e , ip) = 
(F,j3 e ip) =: (F,i/)p), where, for an arbitrary test function -0, ipp is also a test function, since is both 


tempered and C°° (cf. Definition 3.1). Since F is a tempered distribution, this term makes sense. Thus, v, 
is well-defined. 


at 


Unlike in the classical case, the nonlocal antiderivative v e can be much rougher than its “derivative” F. 


This occurs if a t is chosen to be too flat at zero. For example, in the estimate (3.10), if a t is Lipschitz 
continuous at zero (i.e. k a = 1), then we see that *& e (|£|) > C a ^/f; 2 for £ sufficiently large. This means that 
F(£)/6t e (£) ~ F(£)£ 2 for £ sufficiently large, which shows that v e , in this case, is as smooth as F". This gets 
worse the larger that k a is, or the flatter that a e is at zero. 

However, if a e is chosen to be sufficiently singular at zero, then nonlocal integration can be a smoothing 


process, as in the classical case. From (3.10), we see that if k a = — 1 — <5 for 0 < S < 1, then *& e (|£|) > 
Ca,e|£| -<5 - This means that v e ~ F/\£\ S f° r ICI sufficiently large. Therefore, v e is actually smoother than its 
“derivative” F, for these kernels a e . 


Now that we have made sense of the formal solution (3.7), we can address its convergence to the classical 
solution v(t) = f F(t)dt. 


Theorem 3.4 (Weak convergence). Let F satisfy Definition 3-4 and a e Definitions 2.1, 3.1 3.5, and 3.6 
Then v e —v in S 1 as e —>• 0. That is, ( v e ,ip) —> (v,if) for every ip in S. 

Proof. We decompose 2n/id e — 1 /f = jd t into near- and far-held contributions: 

= /5e('S’)l(0,l/2eb c «) (I ^ I) + Pe(^(l/2eb a ,oo) (|£])i 


where b a is as in Dehnition 3.6 Applying the near- and far-held estimates (3.12) and (3.13), we obtain: 
C a e 


&(KI) < 


-l(0,l/2efc o )(|^|) + (^ffi~ e2+ka |£| 1 + fca + 2e6 Q ^ l(l/2e6„,o 


m)- 


2nb a — C a /4b a 

If /3 e represents the nonlocal correction to l/£, then we see that it vanishes pointwise as e —► 0. 


(3.14) 
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We write 2mv e as F /f + Ff3 e and show that the second term vanishes weakly. For any test function ip, 
we let ip := F~ x ip. We have (F/3 e ,ip) = (F,f3 e ip) =: (F,ip e ), where the test function ip e (£) — > 0 uniformly as 
e —> 0. This shows that (v e , ip) = (v e , ip) ( v , ip) as e —> 0. □ 


Remark 3.10. Obviously, we are interested in those kernels a for which the denominator of the first term 


in (3.14) is positive. Such kernels exist. For example, if we assume that supp a is in (—b a ,b a ), then, under 


the hypotheses of Definition 3.6 we can guarantee positivity if —1.402984... < k a < —0.509742.... More 
specifically, if = K + a as well (so that a is a power function), then this holds if —2 < k a < 0 or if 
k a > (24- 27t 2 )/(7t 2 - 6) = 1.1011.... 

If we make a e more singular and restrict A to be a function in an L p space, then we can use the estimate 


(3.14) to obtain stronger convergence results. 


Corollary 3.1 (Strong convergence). Suppose that —2 < k a < — 1, such that a is not integrable at zero. 
Suppose also that F is in L p , where 1 < p < oo. Then v e — v is in L p , and v e — v —> 0 in L p . 


Proof. We have |£| 1+fe “ 
shows that 2iTi(v e — v) 
v £ — v —> 0 in L p . 


< (2MIM- 1 in ( |3.14[ ), which shows that /3 C (|£|) —► 0 uniformly as e —► 0. It also 
= F~ 1 [F'l%] is in L p . Thus, F/3 e —> 0 in L q , where 1/q = 1 — 1/p. This shows that 

□ 
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